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ABSTRACT 

A steradian of the southern sky has been imaged at 151.5 MHz using the Mauritius 
Radio Telescope (MRT) . These images show systematics in positional errors of sources 
when compared to source positions in the Molonglo Reference Catalogue (MRC). We 
have applied two-dimensional homography to correct for systematic positional errors 
in the image domain and thereby avoid re-processing the visibility data. Positions of 
bright (above 15-cr) point sources, common to MRT catalogue and MRC, are used to 
set up an over-determined system to solve for the homography matrix. After correction 
the errors are found to be within 10% of the beamwidth for these bright sources and 
the systematics are eliminated from the images. This technique will be of relevance 
to the new generation radio telescopes where, owing to huge data rates, only images 
after a certain integration would be recorded as opposed to raw visibilities. It is also 
interesting to note how our investigations cued to possible errors in the array geometry. 
The analysis of positional errors of sources showed that MRT images are stretched in 
declination by ~ 1 part in 1000. This translates to a compression of the baseline scale 
in the visibility domain. The array geometry was re-estimated using the astrometry 
principle. The estimates show an error of ^ 1 mm/m, which results in an error of 
about half a wavelength at 150 MHz for a 1 km north-south baseline. The estimates 
also indicate that the east-west arm is inclined by an angle of ^ 40" to the true 
east-west direction. 

Key words: surveys - techniques: image processing - astrometry - techniques: in- 
terferometric - telescope - catalogues 



INTRODUCTION ray (MWA), that are being built are low frequency arrays; 

, , „ , „ , r7r~, ni 1 clearly indicatins a renewed interest in metre-waveleneth 

The Mauritius Radio Telescope (MRT) JColap et al.lll995al : ^, , ^ u • i ■ j ■ ■ i j 

l Udava Shankar et al 1120021 is a Fourier svnthesis T shaned astronomy. The key astrophysical science drivers include ac- 
luuava oiiaiij\.ai ci ai.iizuu A) lo a i^ouiiei byiiLiieoio, _L-Diiapeci i , . , i i i , . • , i i , • • 

; ' celeration, turbulence and propagation m the galactic m- 

non-coplanar array operating at 151.5 MHz. The telescope , . n i. , . ,i i . i i i .n ■ i 

, f ^ , terstellar medium, exploring the high redshiit universe and 

was built to fill the gap m the availability or deep sky sur- j_ ■ j_ ^ n i-r^i i 

,. f . . , , u - i, transient phenomenon, as well as searching tor the red- 



yeys at low radio frequencies in the southern hemisphere 
The aim of the survey with MRT is to contribute to the 



shifted signature of neutral hydrogen from the cosmologi- 

, . , cally important epoch of reionisation (EoR). The surveys 

database oi southern sky sources m the decimation ( ) range , . , .,, . , , . „ . . 

; i , made using such arrays will provide critical mtormation 

—70 to —10 , covering the entire right ascension (a), with , , „ , , .', .„ , . , r i i . i 

, , , about toregrounds which will also provide a useful database 

a synthesised beam of 4 x 4 .6 sec za and an expected point r i . i , i . • i i , ■ ^ T-nm 

_i for both extragalactic and galactic sources. MKi survey at 
source sensitivity (l-cr) of ~ 110 mjy beam . The zenith ^ ,,tt • j_ • j_ j- .l- j ■ i i-^- -n 

, , ^ , ON 151.5 MHz IS a step m that direction and, m addition, wiff 

anqle (za) is given by (o — (p), where, m (« —20.14 ) is . i . i i i i r i-i .■ 

, , . ' r ^ iT^m iT^rr^ , n \ ■ 1 , , providc thc cruciaf sfcy ffiodef for calibratiou. 
the latitude of MKI . MR! has been designed to be the , • . » «-T-.rr^ • .11 1 ,1 ■ 1 

Imaging at MKI is presentfy done onfy on the merid- 

southern-sky equivafent of the Cambridge 6C survey at • ^ • • • 11 r 1 \ ^ 

, ■] „ 7. — . ni K lan to minimise the probfems of non-copfanarity. A two- 

151.5 MHz tealdwin et al.lll985l ). • wo r^^ • • • j- ^ • <■ j 

' \ , T ^ n dimensional (2-D) image m Q-sin2:a coordinates is formed 

The next generation radio telescopes, like the LOw Fre- 1,1. ,• • 1 ■ .1 .1 

.„ „^.„N , , , by stacking one-dimensional (1-D) images on the merid- 

quency ARray (LOFAR) and the Murchison Widefield Ar- . _^ ^ ., w t c ^ j- 

^ ' lan at ditterent sidereal times. Images of ~ a steradian 

(is'' < a 24^^,-70° ^ 5 ^ -10°) of the southern sky, 
with an rms noise in images of ~ 300 mJy beam"' (1- 
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Figure 1. Positional error analysis of ^ 400 sources (above 15-(t) common to MRT catalogue and MRC. For visualisation, the errors 
are shown in percentages of MRT beamwidths. (a) The first row subplot shows errors in a against a; no systematics are observed. The 
second row subplot shows errors in a against sin za; errors show a linear gradient as a function of sin za. (b) The first row subplot shows 
errors in S against a; no systematics are observed. The second row subplot shows errors in 5 against sinza; errors show a linear gradient 
as a function of sin za. The second row subplots in (a) and (b) also indicate declination (or equivalent sin za) ranges of the four zones 
imaged with different delay settings, (c) and (d) show histograms of errors in a and 5, respectively. The histogram of errors in 5 shows 
a broader spread compared to errors in a. 



a), were produced by iPandev and Udava ShankaJ l|2005l l. 
A suite of programs developed in-house was used to reduce 
~ 5000 hours of the survey data (a quarter of the total 
~ 20, 000 hours observed over a span of ~ 5 years). The de- 
convolved images and a source ca talogue of ~ 2, 800 sources 
were published bv lPandevI l|2006l ). 

Systematics in positional errors were found when the 
positions of sources common to MR T catalogue and th e Mo- 
longlo Ref erence C a talogu e (MRC) (|Large et al.lll98ll ) were 
compared. iPandevI (|2006l ) treated the systematics in errors 
in a and sin za independently. By estimating two separate 1- 
D least-squares fits for errors in a and sin za the systematics 
were corrected only in the source catalogue. However, errors 
remained in the images which impede usefulness of MRT 
images for multi-wavelength analysis of sources. In addition, 
the source of errors was not investigated. At MRT, the vis- 
ibility data is processed through several complex stages of 
data reduction specific to the array, especially, arising due 



to its non-coplanarity (|Udava Shankar et al.|[2002l 'l. It was 
therefore decided to correct for errors in the image domain 
and avoid re-processing the visibility data. 

This paper describes the application of 2-D homog- 
raphy, a technique ubiquitous in the computer vision and 
graphics community, to correct the errors in the image do- 
main. Homography is used to estimate a transformation ma- 
trix (which includes rotation, translation and non-isotropic 
scaling) that accounts for positional errors in the linearly 
gridded 2-D images. In our view, this technique will be 
of relevance to the new generation radio telescopes where, 
owing to huge data rates, only images after a certain in- 
tegra tion would be recorded as opposed to raw visibili- 
ties l|Lonsdale et al.ll2009l ). This paper also describes our 
investigations tracing the positional errors to errors in the 
array geometry used for imaging. Our hypothesis on the 
array geometry, its subsequent confirmation endorsed by re- 
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estimation of the array geometry and its effect on the images 
are also described. 

The rest of the paper is organised as follows. Section [2] 
compares positions of sources common to MRT catalogue 
and MRC. The 2-D homography estimation is briefly de- 
scribed in SectionO SectionUpresents the correction scheme 
and typical results. The re-estimation of MRT array geom- 
etry is described in Section (5] Finally, we summarise and 
present our conclusions in Sectional 



2 POSITIONAL ERRORS 

The positions of sources common to MRT catalogue and 
MRC were compared. We used MRC because of its overlap 
with MRT survey, its proximity in frequency compared to 
other reliable catalogues available and, comparable resolu- 
tion (2'. 62 X 2'.86sec(5 -I- 35°. 5)). Moreover, for sources of 
hsted flux density ^ 1.00 Jy (at 408 MHz) the catalogue 
is reported to be substa ntially complete an d, the reliability 
is reported to be 99.9% (|Large et al.lll98ll ). For our further 
discussions, errors in MRC source positions are considered 
random, without any systematics. 

About 400 bright sources common to the two catalogues 
and with flux density at 151.5 MHz greater than 5 Jy (> 15- 
(j) were identified and their positions were compared. The 
sources were labelled as common if they lie within 4' of each 
other. Since MRC has a source density of ~ 0.5 source deg~^, 
the chances of considering two unrelated sources as common 
are extremely low. A flux threshold of 15-(t ensures a source 
population abundant to reliably estimate homography (ex- 
plained in next section). 

The positional errors in a and 5 show no systematics 
as a function of a (refer first rows of Fig. [1^ andfTJa). For 
visualisation, the errors are shown in percentages of MRT 
beamwidths. The errors in a and 5 show a linear gradient as 
a function of sin za. The errors in a, plotted against sin za, 
reach ~ ±10% of the MRT beamwidth (refer second row 
of Fig. [T^). Whereas, the errors in 5, plotted against sin za, 
are significant and reach ~ ±50% of MRT beamwidth. (refer 
second row of Fig. [TJj). Histograms in Fig. [TJ; and Fig. [TJi 
show the distribution of errors in a and 5, respectively. The 
histogram of errors in 5 shows a broader spread compared 
to errors in a. 

Re-imaging, to correct for errors in the images, would 
involve re- reducing the ~ 5, 000 hours of observed data. Ow- 
ing to the complexity involved it was decided to correct 
for the positional errors in the images, thus avoiding re- 
processing. The 2-D homography estimation technique was 
employed for correcting positional errors in images and is 
discussed in detail in the following section. 



3 2-D HOMOGRAPHY 

The 2-D planar homography is a non-singular linear re- 
lationship between points on planes. Given two sets of 
K corresponding image points in projective coordinates, 
(pt, an d Pfc) € P^, homography map s p^, to the correspond- 
ing pj. l|Hartlev and ZissermanlboOOl 'l . Where, k = 1, . . . , K . 
The homography sought here is a non-singular 3x3 matrix 



H such that: 
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Where, {Xf.,y^) and {x'f.,y'^) represent (a,sin2:a) of K cor- 
responding MRT and MRC sources, respectively. 

In Equation [TJ (ij,, j/j,, 1) and {x'^,y'^, 1) are referred to 
as the homogeneous coordinates and are always represented 
one dimension higher than the dimension of the problem 
space. This is a commonly used representation in computer 
graphics. The simple reason is that with a 2 x 2 matrix 
one can only rotate a set of 2-D points around the origin 
and scale them towards or away from the origin. A 2 x 2 
matrix is incapable of translating a set of 2-D points. The 
homogeneous coordinates allow one to express a translation 
as a multiplication. A single 3x3 matrix, with homogeneous 
coordinates, can account for rotation, scaling and translation 
of 2-D coordinates. For example, from Equation [H a;*, = 
hiiXi, + hi2yk ± ^13- Notice, hi3 (representing translation 
in a-dimension) is simply being added to the normal dot 
product (/iiiXfe ± hi2yk) that together represents rotation 
and scaling. In homogeneous coordinates, the 2-D problem 
space is a plane hovering in the third dimension at a unit 
distance. 

A general homography matrix, for projective transfor- 
mation, has 8 degrees-of-freedom (DOE). For our system, 
both errors in a and 6 have only sin 2;a-dependency. There- 
fore, a less general, 2-D afiine transformation is sufficient. 
A 2-D afhne transformation (two rot ations, two translations 
and two scalings) requires 6-DOF (|Hartlev and ZissermanI 
l2000l ). therefore in H, /131 and /132 are zero. Since each 2-D 
point provides two independent equations, a minimum of 3 
point correspondences are necessary to constrain H in the 
afiine space. A set of K such equation pairs, contributed by 
K point correspondences, form an over-determined linear 
system: 
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h= [ hii,hi2,hi3,h2i,h22,h23,h3i,h32 ]^and. 



b = 
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) XkjVk 
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In Equation [21 T represents transpose of a matrix. This sys- 
tem can be solved by least squares-based estimators. 

At this stage it is useful to consider the effect of us- 
ing (a, sin 2;a)-coordinates to represent the brightness distri- 
bution on the celestial sphere. Ideally, it is the directional 
cosines {I, m, n), with respect to the coordinates of the array, 
which represent the spherical coordinates in the sky. There- 
fore, the image coordinates in which homography should in 
principle be estimated are (Z,m). However, at MRT, for 1-D 
imaging on the meridian: 

m = sin za. (3) 
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Figure 2. Block schematic of tlie correction scheme. Rectangular boxes represent processes; rounded boxes, data or results. 



Therefore, sin za is a natural choice for one of the coordi- 
nates and is indeed used in the present case. On the merid- 
ian, the directional cosine I is zero. For small errors, Al, in 
I, i.e. close to the meridian: 



Al — cos 5 Aa. 



, Aa — Al sec 5. 



(4a) 



(4b) 



Here, Aa is the error in a. Equation |4] shows that an error 
in I will lead to an error in a with a sec 5-dependence. 

The 2-D images of MRT are 1-D images on the merid- 
ian made at different sidereal times and stacked. There- 
fore, positional errors both in a and S do not show sys- 
tematics as a function of a {first rows of Figs. [T^ andfTJa). 
We preferred (a, sin 2a)-representation because all MRT im- 
ages were already generated in this coordinate system. This 
choice compelled us to seek solutions for errors in a as a 
function of sin za rather than sec 5. We plotted errors in a 
against both sec S and sin za and obtained separate linear 
least-squares fits. The rms of residuals in both fits is ~ 5% 
of the beamwidth in a. However, the rms of difference be- 
tween the fitting functions sec S and sin za in the 5 range of 
MRT (-70° to - 10°) is only ~ 1.5% of the beamwidth in 
Q. Therefore, the random errors in the source positions are 
larger than the errors introduced by the preferred (a, sin za)- 
coordinates for pj. and pj.. 

In p^. and p^, the a ranges from 18 hours to 24 
hours and the sin za ranges from —0.8 to 0.2 (correspond- 
ing to the declination range of —70° to —10°). Moreover, 
in matrix A (refer Equation [2]) there are entries of I's 
& O's. Such a matrix is ill-conditioned and in the pres- 
ence of noise in the source positions, the solution for an 
over-det ermined system may diverge from the correct es- 
timate (|Hartlev and Zissermanlbood ). The effect of an ill- 
conditioned matrix is that it amplifies the divergence. A 
normalisation (or pre-conditioning) is therefore required. 



3.1 Data normalisation and denormalisation 

To obtain a good estimate of the transformation matrix 
we ad opted the normalisation scheme proposed bv iHartlevI 
The normalisation ensures freedom on arbitrary 
choices of scale and coordinate origin, leading to algebraic 
minimisation in a fixed canonical frame. The homography 
matrix H is estimated from normalised coordinates by the 
least-squares method using singular value decomposition 
(SVD). The matrix is then denormalised to obtain H. The 
scheme is briefly described below: 



(i) Normalisation of p: Compute a transformation ma- 
trix M, consisting of a translation and scaling, that takes 
points pj. to a new set of points pj. such that the centroid 
of the points pj. is the coordinate origin (0,0)"^, and their 
average distance from the origin is V2. 

(ii) Normalisation of p': Compute a similar transfor- 
mation matrix M', transforming points pj. to p'^.. 

(iii) Estimate homography: Estimate the homography 
matrix H from the normalised correspondences f>^. p'^. 
using the algorithm described earlier in the main section. 

(iv) Denormalisation: The final homography matrix is 
given by: 



H = M 



HM. 



4 THE CORRECTION SCHEME 

Fig.[2]shows the block schematic of the correction scheme. At 
MRT, the full declination range for each sidereal hour range 
is divided into 4 zones (refer second row in Fig. [1^ orfTJj). 
Each zone is imaged with different delay settings to keep the 
bandwidth decorrelation to < 20%. Therefore, the 6 side- 
real hours of images under consideration, have 24 images 
(~ 15° X 15°). 

Using the population of common sources, there are four 
possible alternatives to correct MRT images by computing: 

(i) 24 homography matrices - one for each image. 

(ii) 6 matrices - one for each sidereal hour. 

(iii) 4 matrices - one for each declination zone. 

(iv) A single homography matrix for the entire steradian. 

In principle, bright sources in each image (15° x 15°) can 
be used to independently estimate a homography matrix. 
Our earlier experiments to correct each image independently 
showed that the homography matrices were similar. The 
plots of errors in a and S plotted against a and sin za (refer 
to Fig. [1^ and [T]d) indicate that the errors are independent 
of the four delay zones and the range of a. This implies that 
estimating a single homography matrix for the entire source 
population should suffice in representing the errors. 

The homography matrix estimated using ~ 400 com- 
mon sources (described in Section [2} is: 



H = 



1.0000 0.0006 
0.0000 0.9990 
0.0000 0.0000 



0.0001 
-0.0009 
1.0000 



(5) 



In the estimated homography matrix, hn = 1.0000 indicates 
there is no correction required in a as a function of a. hi2 = 
0.0006 indicates MRT images should be corrected in a with a 
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Figure 3. Positional error analysis after homography-based correction, (a) First and second row subplots show errors in a against a 
and sin za, respectively, (b) The first and second row subplots show errors in 5 against a and sin za, respectively, (c) and (d) show 
histograms of errors in a and &, respectively. A comparison of these plots with Fig. [T] demonstrate that homography has removed the 
systematics and the residual errors are within 10% of the beamwidth. 
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Figure 4. Scatter plot of errors in a and <5. (a) Before correction and (b) after homography-based correction. After correction the scatter 
is almost circular as opposed to elliptical before correction. The rms before correction is 20% of the beamwidth. After correction, the 
rms is reduced to 7% of the beamwidth and, the systematic errors have been removed. 
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Figure 5. MRT contours overlaid on SUMSS image, (a) and (c) MRT contours before correction for sources at about S 



-66° 



and 5 = —40°, respectively, (b) and (d) Corresponding MRT contours after homography-based correction show that 2-D homography 
corrected the positional errors. Notice, (c) and (d) are included here for visual emphasis. Since the errors around S = —40° are within 
10% of the beamwidth the contours show a good overlap both before and after and, as expected homography has not applied perceivable 
correction to images at this declination. 



sin za dependence. The estimated correction is up to ~ 10% 
of the beam in a, at the extreme ends of the sinjza range. 
Similarly, /121 ~ 0.0000 indicates that there is no correction 
required in sin za, as a function of a. However, /122 ~ 0.9990 
indicates that MRT images should be compressed in sin za 
by a factor of 0.9990 (which is ~ 1 part in 1000). The values 
of and /i23 indicate that the zero cross-overs of errors in 
both a and sin za plotted against sin za are close to the sin za 
of the calibration source (MRC 1932-464) used for imaging. 

Using Equation [T] the homography matrix is used to 
project each pixel from the images to a new position, effec- 
tively correcting for positional errors in images. 



4.1 Corrected images and discussion 

Fig. [3] shows positional errors in 5 after homography-based 
correction. A comparison of these plots with Fig. [1] demon- 
strate that homography has removed the systematics and 
the residual errors are within 10% of the beamwidth for 
sources above 15-a, as expected. Fig. [4^ and Fig. |4|3 show 
scatter plots of errors in S against errors in a before and af- 
ter correction, respectively. For visualisation, the errors are 
represented in percentages of respective MRT beamwidths. 
Notice, after correction (refer Fig. the scatter is al- 



most circular as opposed to elliptical before correction (re- 
fer Fig. [4^). The rms before correction is ~ 20% of the 
beamwidth. After correction, the rms is reduced to ~ 7% 
of the beamwidth and, the systematic errors have been re- 
moved. 

Fig. [5^ and [Sja show MRT contours before and after 
correction, respectively, overlaid on S UMSS (Sydney Un i- 
versity Molonglo Sky Survey) image (|Mauch et all I2003I '). 
for a source around S = —67°. The corrected MRT image 
contours in Fig. \Sjp overlap with the source in SUMSS im- 
age. Figs. [5j; and [5]i show similar comparison for a source 
around S — —40°. Notice Fig. [TJi, since the errors around 
5 — —40° are within 10% of the beamwidth, the contours 
in both Figs. \5jf and (5)1 show a good overlap as expected 
and homography has not applied perceivable correction to 
images at this declination. We have overlaid MRT contours 
on a number of extended sources at 843 MHz reported by 
Ijones and McAdamI (|l992l ) . Fig. [6] shows a typical overlay 
of MRT contours on SUMSS image of a region around the 
cluster Abell 3667. The overlay is perceivably satisfactory. 

The 2-D homography corrected the positional errors in 
the image domain. For imaging the remaining ~ 3.5 steradi- 
ans of MRT survey, ~ 15000 hours of data has to be reduced. 
Ideally, for imaging the new regions, one would like to trace 
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Figure 7. Typical calibrator phase differences (in radians) of 
MRC 0915-118 & MRC 1932-464, plotted against v (in metres). 
The straight line shown is a linear robust fit obtained for the 
data. 



the source of these errors and correct them in the visibihties. 
In the following section we discuss how we traced the source 
of errors and corrected them in the visibility domain. 



5 ARRAY GEOMETRY: HYPOTHESIS & 
RE-ESTIMATION 

This section describes our expansion- compression hypothe- 
sis for the source of errors in our images. The subsequent 
corrections we estimated and applied to eliminate the errors 
are also described. 

For meridian transit imaging, m = sin za. The bright- 
ness distribution in the sky as a function of sin za and 
the complex visibilities measured for different values of 
the north-so uth (NS) baseline vector compo nent v form a 
Fourier pair (|Christiansen and Hogbomll 19851 ). A scaling er- 
ror of K in m will result in a scaling factor of in the 
?;-component of the baseline vector. By positional error anal- 
ysis it is clear that MRT images are stretched {expanded) in 
declination, i.e.. 



^cd ^ "^true 



, V. 



measured 



^true ■ 



(6a) 



(6b) 



Note, for images the 2-D homography estimated a cor- 
rection {compression) factor, , of 0.9990. This cued to 
the hypothesis that we have compressed the north-south 
baseline vectors. Equation [Sja means, a baseline distance of 
~ 1000 m in the NS arm was wrongly measured as ~ 999 m 
(1 part in 1000). Similarly, a sin za-dependent correction in a 
cued to possible ii-component in the east-west (EW) base- 
line vectors. Next, we describe the re-estimation of array 
geometry. 

We begin with a brief description of the mode of obser- 
vations with MRT. MRT has 32 fixed antennas in the EW 
arm and 15 movable antenna trolleys in the NS arm. For 
measuring visibilities, the 15 NS trolleys are configured by 
spreading them over 84 m with an inter-trolley spacing of 



6 m (to avoid shadowing of one trolley by another). MRT 
measures different Fourier components of the brightness dis- 
tribution of the sky in 63 different configurations (referred to 
as allocations) to sample NS baselines every 1 m. Therefore, 
effectively, there are 945 antenna positions (63 allocations * 
15 antennas/allocation) in the NS arm and a total of 30,240 
(945 * 32) visibilities are used for imaging. 

A small error in a measuring scale of relatively shorter 
length is likely to build up systematically while establishing 
the geometry of longer baselines. This effect would be ob- 
served in the instrumental phases estimated using different 
calibrators. In principle, the instrumental phases estimated 
using two calibrators at different declinations, for a given 
baseline, should be the same, allowing for temporal vari- 
ations in the instrumental gains. A non-zero difi'erence in 
these estimates may be due to positional errors of the base- 
line or positions of calibrators. As mentioned earlier, our 
analysis of positional error in sources and the homography 
matrix cued to positional errors in baselines (or antenna po- 
sition s) . The simple principle of astrometry ({Thomson et al.l 
[2003) was used to estimate errors in antenna positions and 
is discussed below. 

The observed visibility phase, tp^^ , in a baseline with 
components {uij,Vij,Wij), due to calibrator Si with direc- 



tion cosines {r^ , m 



V'j^.^ = l^^ Uij + rri 



^), is given by: 
•Si 



<ms 



(7) 



Where, (j)™ represents true instrumental phases, i = 
1, 2, . . . , 32 represents EW antennas and j — 1,2,..., 945 
represents NS antennas. For meridian transit imaging Equa- 
tion [7] becomes: 



= -v,j sin (^Z^'^i j + cos (^ZA^^^^ + 4>. 



ms 
ij ■ 



(8) 



The instrumental phases, (j)^^ , estimated using the measured 
geometry are given by: 

(?!)f/ = -A«,j sin [ZA^^^ + Awij cos (^ZA^^^ + ^jf . (9) 

Here, Aw^ and Aw^j are errors in the assumed baseline 
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Figure 8. (a) Estimate of the antenna position along the North-South arm of MKT array. The fit to the estimates shows a gradient of 
1 part in 1000, along the north-south, (b) Estimate of the error in v coordinate of MRT east-west arm. The fit to the estimate shows a 
gradient of 2 part in 10,000, along the east-west. 



are phases of complex baseline gains obtained 



vectors 

in the process of calibration. Equation O has three un 
knowns. To reduce the number of unknowns, one can elim- 
inate the true instrumental phases by taking a difference 
= 'Pij ~ '/'i/ ) between the instrumental phases es- 
timated using two calibrators. This difference gives: 



; l^sin ^ 
l^cos ^ 



ZA^ij - sin [^ZA^ 
ZA^A - cos (ZA^ 



')] 



(10) 



Note, the ui-components of the baseline vectors are short 
and non-cumulative measurements. Therefore, in principle, 
one can consider Aw^j as zero-mean random errors with no 
systematics. Equation [TU] in that case can be written as: 



A<;/)f/^ = -Av,^ [sin (^ZA^^'^ - sin (^ZA^ 



(11) 



Describing the system in terms of errors in antenna posi- 
tions, as opposed to errors in baseline positions, Equation llll 
becomes: 

,5, 



AC 



{Av, - Avj) [sin ( ZA 



sin ZA" 



(12) 

This equation is also not sufficient to solve for errors in the 
antenna positions as we have two unknowns and one equa- 
tion. We set up another equation using a third calibrator 
source, 53, spaced away in declination from Si and 52: 



Av, - Avj) [sin ( ZA^A - sin ( ZA^^ 



(13) 

The Equations [12] and [13] are a linear set of equations for 
one baseline. For the measurements in 63 allocations, the set 
of equations can be formulated in a matrix form and solved 
by SVD-based least-squares estimator: 

Avi = b. (14) 

Where, the measurement vector x £ R'^ is to be determined. 
Here, c = 977. The measurement vector gives Avi and Avj 
estimates for 32 EW and 945 NS antenna locations, respec- 
tively. The observation vector b consists of two sub- matrices, 



bi £ R' 1 and b2 G R'^ , formed using the left-hand-side of 
Equations [12] and 1131 respectively. Here, r^ = r2 = 30240, 
i.e., the total number of visibilities measured for imaging. 
Therefore, b e R'^^^^o. The data matrix A £ r60480x977^ 
Each row in the data matrix has only two non-zero elements, 
corresponding to a baseline formed by one EW and one NS 
antenna, making it very sparse. 

The observation vector is constructed from the gain ta- 
bles of the array obtained using calibrators MRC 0407-658 
(5i), MRC 0915-118 (52) and MRC 1932-464 (^3). The sen- 
sitivity per baseline at MRT is ~ 26 Jy for a 1 MHz band- 
width and an integration time of one second. It takes ~ 10 
minutes of time for sources at 5 = —40° to transit a 2° 
primary beamwidth of elements in the east-west array. This 
leads to a sensitivity per baseline (including the non-uniform 
weighting due to primary beam) of ~ 2 Jy. The flux density 
of these three calibrators as seen by MRT is ~ 100 Jy; strong 
to get reliable calibration. Further, the calibrators are un- 
resolved and isolated from confusing so urces and have well 
known measured positions l|Goladl 19981 ). 

A plot of typical phase differences obtained using the 
pair of calibrators S2 and ^3 is shown in Fig. [7] Fig.[8K shows 
the estimated errors in 945 NS antenna positions. The errors 
show a gradient of 1 part in 1000 along the NS arm. This 
matches with the linear gradients in the phase differences es- 
timated from the calibrators. The estimates in Fig. [8]3 show 
alignment errors of the 32 antennas in the EW arm along 
the NS-direction. The fit shows a gradient of about 2 part in 
10,000. This indicates that the EW arm is mis-aligned from 
the true EW-direction. At one extreme end (1 km from the 
centre of the array) of the EW arm the error is ~ 0.2 m, 
equivalent to an angular distance of ~ 40" from the centre 
of the array. This is the source of a small sin za-dependent 
error in a that was observed in both positional error analy- 
sis and the homography matrix. Further, our simulation of 
the synthesised beam in a with old EW antenna positions 
and the corrected EW antenna positions indeed confirm this 
sin 2;a-dependent error in a. 

Using the new antenna positions we have re-imaged one 



2-D homography-based correction of positional errors in widefield MRT images 9 



hour from the steradian and have also imaged a completely 
new steradian. We find no systematics in positional errors 
thus endorsing our re-estimated array geometry. 



6 CONCLUSIONS 

The homography-based correction was able to correct for 
systematics in positional errors in the image domain and the 
errors are within 10% of the beamwidth for sources above 
15-cr. The corrected images of one steradian are available for 
download at http://www.rri.res.in/ surveys / MRT\ 

Positional error analysis showed that uncorrected MRT 
images are stretched in declination by ~ 1 part in 1000. 
This translates to a compression of the NS baseline vec- 
tor, in the visibility domain. The analysis also showed a 
sin za-dependent error in a. This cued towards possible er- 
rors in our estimation of the array geometry. By formulat- 
ing a linear system, using instrumental phases estimated 
from three well separated calibrators whose positions are 
well known, the array geometry was re-estimated. The esti- 
mated error in the u-component of the NS baseline vectors 
is about 1 mm/m. In other words, the error is about half 
a wavelength at 150 MHz (1 m) for a 1 km baseline. The 
estimates also show a small (2 part in 10,000) ii-component 
in the purely EW baseline vectors. This indicates that the 
EW arm is mis-aligned and inclined at an angle of ~ 40", 
to the true EW direction. These estimates match with the 
observed stretching of MRT images shown by both the po- 
sitional error analysis and the homography matrix. 

Using the new antenna positions we have re-imaged 
one hour from the steradian and have also imaged a com- 
pletely new steradian. We find no systematics in positional 
errors. This endorses our re-estimated array geometry. Re- 
imaging one steradian starting from visibilities would have 
been a very time consuming exercise. Development of 2-D 
homography-based correction enabled us to correct for the 
positional errors in the image domain. In our view, this new 
technique will be of relevance to the new generation radio 
telescopes where, owing to huge data rates, only images af- 
ter a certain integration would be recorded as opposed to 
raw visibilities. 
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